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Abstract 

The Kastor-Traschen metric is a time-dependent solution of the Einstein-Maxwell equations 
with positive cosmological constant A which can be used to describe an arbitrary number of 
charged dynamical black holes. In this paper, we consider the null geodesic structure of this 
solution, in particular, focusing on the projection to the space of orbits of the timelike conformal 
retraction. It is found that these projected light rays arise as integral curves of a system of 
third order ordinary differential equations. This system is not uniquely defined, however, and 
we use the inherent freedom to construct a new system whose integral curves coincide with 
the projection of distinguished null curves of Kastor-Traschen arising from a magnetic flow. 
We discuss our results in the one-centre case and demonstrate a link to conformal circles in 
the limit A — > 0. We also show how to construct analytic expressions for the projected null 
geodesies of this metric by exploiting a well-known difieomorphism between the K-T metric and 
extremal Reissner-Nordstrom deSitter. We make some remarks about the two-centre solution 
and demonstrate a link with the one-centre case. 



1 Introduction 

Recently, the study of the behaviour of null geodesies has led to the development of geometric tools 
useful for revealing physical properties of spacetimes which admit a particular type of timelike 
vector field. 

In the case of a static metric, for example, a convenient geometric tool is the optical metric defined 
on the space of orbits of a given hypersurface-orthogonal (HSO) Killing vector. Light rays of 
the static metric are found to project down to unparametrised geodesies of the associated optical 
metric. This property has been used to understand some observations of the dynamics of light rays 
in Schwarzschild-deSitter metrics [T] , and to give an alternative interpretation of black- hole no-hair 
theorems [2]. Static metrics admitting more than one HSO Killing vector have been considered in 
the context of projective equivalence of the associated optical metrics [3]. 

Some progress has also been made in the case of stationary metrics [4], where there are two 
distinguished geometric structures on the space of orbits B of the timelike Killing vector field. 
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On one hand, the light rays of the stationary metric project down to solutions of the Zermelo 
navigation problem on B with specially-defined metric h and wind W. Alternatively, one can 
construct a Randers-Finsler structure on B with Randers data {ay, h} such that the null geodesic 
flow projects down to the magnetic flow due to db on the curved manifold {B, Cbij}. 
In this paper, we generalise this work to another class of metrics - those which admit a timelike 
conformal retraction, i.e, there exists a hypersurface-orthogonal timelike vector field, O, for which 
the conformal structure on its space of orbits is preserved along the integral curves of 9. In the 
Riemannian case, such metrics have arisen as supersymmetric solutions of minimal Af = 2 gauged 
supergravity with anti-self-dual Maxwell field [5]. Moreover, when the anti-self-duality condition 
is relaxed in the case of positive cosmological constant, one obtains a solution, also admitting a 
conformal retraction, which is the Riemannian analogue of the well-known Kastor-Traschen metric 

In Lorentzian signature, this metric is a time-dependent solution to the Einstein-Maxwell equations 
which can be seen to describe an arbitrary number of dynamical charged black holes in a deSitter 
background. In this paper, we focus on the null geodesic structure of these metrics and in particular, 
it is found that light rays project down to integral curves of a system of third order ordinary 
differential equations (ODEs) on the space of orbits of 0. In particular, we find that the projected 
null geodesies form only a subset of the set of integral curves and that this yields a freedom in how 
this third order system is defined. 

In section 4, we pay particular attention to the one-centre solution (single black hole). We show 
that, in the limit as the cosmological constant tends to zero, our system of ODEs becomes that 
describing conformal circles of the flat metric (as described in [7]). We use this result to motivate 
section 5 where we discuss the numerology of the problem and construct a new third order system 
in three dimensions. The advantage of this new system is that it will allow us to give a physical 
interpretation to those integral curves which do not arise as projected null geodesies, coinciding 
with motion in the background magnetic field. This formulation also allows us to characterize the 
projected light rays. Proposition 5.1 is the central result of this paper. 

In section 6, we show how to derive analytic expressions for projected null geodesies of the one- 
centre case by making use of a diffcomorphism between the extreme Reissner-Nordstrom deSitter 
metric and the Kastor-Traschen metric. We give plots of some of these curves and discuss the 
horizon structure in both sets of coordinates. 

In section 7, we discuss the Kastor-Traschen solution with two centres. We demonstrate that, for 
special initial conditions, the projection along the conformal retraction of a null geodesic will lie in 
a plane. For one such curve, we illustrate a connection between the null geodesies of the two-centre 
Kastor-Traschen metric and a third order system that arises in the analysis of the one-centre case. 
We also look at the perturbations away from this plane and give a strict condition for stability. 
Throughout this paper, we refer to the projection of geodesies onto the space of orbits of the 
conformal retraction O as the retraction projection to avoid ambiguity in later sections. 
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2 Conformal Retraction in the Kastor-Traschen metric 



The Kastor-Traschen metrics are a class of time dependent solutions to the Einstein-Maxwell equa- 
tions with positive cosmological constant A [B]. In local coordinates, these metrics may be written 

as 

dT 2 



(v + cry 

with Maxwell 1-form 



ff=- ^T^ 2 +{V + cTfh (1) 
dT 



V + cT' 

where V = V(x) is a harmonic function on the spatial coordinates, c = ±y^ is a constant and 

h = hijdx l dx 1 , i,j = 1, 2, 3 

is the flat 3-dimensional Riemannian metric (which we express in Cartesian coordinates for now). 
Notice here that since V is a harmonic function there is a freedom in defining the electromagnetic 
field tensor F = dA so that the Einstein-Maxwell equations are still satisfied i.e, we can write 

F = L= (dA - ue ijk h il ^ i dx j A dx k ^j , i,j,k= 1,2,3 (2) 

where v is a constant and e^-fc is totally anti-symmetric in its indices with €123 = 1. This allows us 
to introduce a magnetic field B oc W into our definition of the Kastor-Traschen solution, a notion 
which will be useful later. 

In the limit as c — > 0, these metrics reduce to the well-known Majumdar-Papapetrou solutions jS], 
[5]. For the M-P metrics, it can be shown that solutions with 



X - W n 



where |x — w Q | = (hjk(x? ' — w 3 a )(x h — w^)) 1 / 2 and w a is a fixed vector for each a, can be analytically 
extended to be interpreted as a system of charge equal mass black holes [TU]. In this system, the 
gravitational forces on each black hole are balanced by the electrostatic forces. However, when 
c 7^ 0, the black holes are dynamic and can be observed to coalesce [6j. 

As noted in the introduction, an interesting property of the K-T metrics is that they admit a 
timelike conformal retraction 0. In other words, 

£-eHnv = fHnv + @(^Ci/j 

where 

and / and C are an arbitrary function and one-form, respectively. Greek indices here run over 
the values 0,1,2,3 and we raise and lower indices using the metric g^ v . In our coordinates Q, the 
conformal retraction is = Furthermore, the tensor H is given by 

H 0fl = , H^y = g^ y otherwise. 
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The Lie derivative of this tensor is easy to compute 
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(C e H)^ = Q X H^ X + H^O^ + H^e^ = —H^ = 2c(V + cT)H^. 

Hence, in this case, our function / and one-form C are given by 

f = 2c{V + cT) , C = 0. 



Now, let us decree that under a change of metric g = Org, the choice of conformal retraction 

Then, H„ v = tt 2 H^ 



remains unchanged i.e, = = Sp. Then, = fl 2 !!^ and 



(£e#V = {Ce0 2 H), lv = (£ e n 2 )H^ + n 2 (£ e H)^ = (/ + 2£ e logO)^ + Q 2 Q (fJl C u) . 
Hence, 

/=/ + 2£ logn , C = tl 2 C. 
and, for our example, C = for any choice of metric in the conformal class of g. 

3 Projection of Null geodesies with arc-length parametrisa- 
tion 

Since we are considering only the null geodesic structure of the Kastor-Traschen metrics, we may 
as well begin with a conformally rescaled version of - this will ease the computation a little. So, 
let us take a new definition of g 

1 dT 2 ■ k 

9 {V + cT) 29 = -{V + cTY + h * dX dX ■ 

For this metric, we can calculate the Christoffel symbols 



■ dV 

r« — M.7 r* — n — v l 



Hence, the geodesic equations for the spatial components may be written as 



(V + cT) 5 dx l w > ds 

where F(s) is some function of the curve parameter s. If we assume that our geodesies are null, 
then we can rewrite this equation as 

h a — l h jk x^x k = F{s)x\ (3) 



V + cT dx 

Now let us impose the condition that s be the arc-length parameter for the metric h, i.e, hj^x 3 x k 
1. 

w V + cTdx k 
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Hence, given that we use the arc-length parametrisation, the equation for null geodesies ^ may 
be written as 



V + cT V dx l dx k 
In the case when h in ([T]) is non-flat, we can derive an analogous expression 



2 ay 

where r* fc are the connection components of the metric /i. Two properties that we get from equation 
Q are 

fci** 7 '** = tFTvtS** ( 6 ) 



y + cT <9x fe * 
<9x fc 7 + I \dx k )\dx k ) \dx k 



(7) 



Now, if we differentiate Q, we get the following 



2 , A f u dV dV . k ., 

-x + cl h — -j — — — j- x x 



(V + cT) 2 \dx k J V dx l dx k 

2 ( h u d 2 V ±m _ d 2 V ±k . l±i _ d^ ±k±i _ dV^ 
V + cT \ dx l dx m dx k dx l dx k dx k 



Then, using the null condition to eliminate T and equations (|4j), (|6| and (|7j) to eliminate the T 
parameter, we can rewrite this in the form: 



3x fc !^ dV , ( ,dV dV t 

., T-> r --T- — — ir l T l - 2c h %l — r k r' 

( 9V dV _ (9V . k \2\ dx l \ dxl dx k x x 

\dx k dxk \dx k / J 



2 

2 ( d 2 v .„ n-v 



V + cT V" dx l dx m ~ dx k dx lXXX )' ^ 

Here, we can also eliminate the factor of v + cT by using (7| so that we really have a system of third 
order ODEs completely dependent on the spatial coordrnatcs of g only. However, it is useful to 
keep it in the form above for the computation in the next section. 

4 One-Centre Case 

As was described in [5J, a spacetime containing N charged black holes with masses m a (a = 
1, . . . , N) and charges q a = m a in a deSitter background can be represented by equation with 

N 



t-^t V — w 



a=l 1 c 



where w Q is a fixed vector for each a. It is easily verified that V is a harmonic function. 

In this section, we look at the case N = 1 where the black hole is situated at the origin. So, put 



5 



V — r^r. With this definition, we can obtain the following identities: 

|x| ' o 

dV m i d 2 V _ 3?Ti m m 

9^ ~ ~\^ hjlX ' d^fe* ~ jxf^ kmX ~ W jk ' 

In particular, equation Q becomes 

where x.x = hjkX^x k and all subsequent dot products are taken with respect to the metric h unless 
otherwise stated. This allows us to reduce the last term of ffify, that is 



V + cT \ dx l dx m dx k dx l 

2 /3m i 3m . 2 .A 3 { 

-(x.x)x _ — (x.x) x = - — ( x .x)a: 



y + cT vi x l 

With this simplification in mind, we can rewrite our system of third order ODEs ^ as 

x< . -HHV + ^ - (x.x)*') - 3(x.x) (Jp + 2(|x|i ^ xx)2) ) *< (10) 

Any null geodesic of the Kastor-Traschen metric g will project down to an integral curve of this 
system of third-order ODEs. 

4.1 Conformal Circles 



As c — > 0, it's obvious that the second term on the right-hand side of ( 10 ) vanishes. However, using 
the second equation in Q with c = and V = j^, we also find that 

x.x=2i&s:-i 



and the third term vanishes. To see the vanishing of the third term explicitly occuring with the 
vanishing of c, it seems we need to reintroduce the time coordinate T, in some way. For example, 
using ([7| we can write ( 10 1 as 

x* = -|x| 2 i l + ^ {x l - (x.x)i l ) + 3(x.x) ( cT l x K x - x ) 



V2m(|x| 2 - (x.x) 2 ) 



Hence, as c — > 0, null geodesies satisfy 

x l = -|x| 2 ±\ (11) 



We shall see that this system (111 occupies a central role in the theory of conformally flat manifolds. 
In general, given a conformal structure [h] on an n-dimensional manifold, there is a distinguished 
family of curves, known as the conformal circles. These curves arise as the integral curves of a 
system of third order ODEs, see [11]. To write this system down, let us choose a metric h in the 
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conformal class with torsion-free connection T 
Then, the Schouten tensor P.\ k is defined as 



FJjfc)) Ricci tensor Rj k and scalar curvature R. 



x i +T i jk x j x k so that U.A = h jk WA k , 



Furthermore, define the vector components U l = x 1 and A 1 

etc, and scalar products coincide with our previous definitions when h is flat. Then, a curve is a 
conformal circle of [h] if it satisfies 



ds 



3|Af 
2IUI- 



(U^tPPj - 2P jk U J U k U l 



(12) 



where there is no restriction on the parameter s. Equation (12) is invariant with respect to confor- 
mal transformations, h — > £l 2 h, and so, conformal circles are defined invariantly by any metric in 
the conformal class. These curves have arisen in a physical context in |12) where the authors have 
used them to discuss the asymptotics of Einstein's equations. Furthermore, properties of "confor- 
mal geodesies" (lifts of conformal circles to the bundle of second order frames over the manifold 
endowed with the conformal Cartan connection) in vacuum and warped-product spacetimes have 
been studied in [13] . 

It is shown in [11] that the conformal circles of a given conformal manifold can be equally defined 
as the set of integral curves of the system of ODEs 



ds 



t} k A j U k = -|A| 2 t7* 



WP* - P jk U J U k U 



(13) 



where, here, s is required to be the arc-length parameter of the metric h. This formulation was 
originally given by Yano in [7] and is more useful for our purposes. If we now let h be the flat 
Riemannian metric, equation (13) reduces to (11). 

Hence, as c — > in the one-centre Kastor-Traschen metric, null geodesies project down to conformal 
circles of the flat metric in three dimensions. It is easily verified, that the set of integral curves of 
(11) in three dimensions is precisely the set of all circles in K 3 . 

We should note here that we get the same result when we let m — > 0. Again, the second term on 
the right-hand side of ( 10 1 obviously vanishes and ^ reduces to 



x l = 0, 



so the third term also vanishes. 



5 Characterisation of Null geodesies 

In the preceding sections, we have determined a system of third-order ODEs which the projected 
null geodesies of g along the conformal retraction 9 must satisfy. However, it is not clear that, 
given an integral curve of Q, it will necessarily be the projection of some light ray of g. In fact, 
we can show this not to be the case and it transpires that the third order ODE system |8]) is not 
uniquely defined. In this section, we discuss this point and construct a new third order system in 
three dimensions for which the integral curves constitute a retraction projection of a special set of 
null curves of the Kastor-Traschen metric, which have a physical interpretation. The projected null 
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geodesies form a subset of these curves which we can characterize. 



For example, let us consider the case c —¥ for the one-centre metric. Here, we found that 
the integral curves of ( 11 1 will be the set of circles in IR 3 . 



However, as c — ► 0, the metric g becomes static with static Killing vector 8 = ^ and it is well 
known that the null geodesies of this metric project down to the unparametrised geodesies of the 
associated optical metric 

/ to \ 4 
h op t = f rr J h jk dx j dx k . 

One can check that the unparametrised geodesies of this metric will be precisely the set of circles in 



R 3 which pass through the origin. Hence, only a subset of the integral curves of ( 11 ) will coincide 
with the projected null geodesies of g. Note that this is also the case for to — > where projected 
null geodesies are described by x 1 = - straight lines. We can check the general numerology here 
to see what happens. 



Firstly, we note that the set of unparametrised geodesies of an arbitrary metric g on some open set 
U C E™ will lift to a foliation, by the geodesic spray, of the projectivised tangent bundle F(TU) 
which we can think of as a 1-dimensional fibration over some (2n — 2)-dimensional space, Z, with 
each point in Z coinciding with a unique geodesic in U. Hence, the set of unparametrised geodesies 
of an n-dimensional manifold constitute a (2n — 2)-parameter family of curves. 
Taking the specific example of the Kastor-Traschen metric g, we have n — 4 and so the number 
of parameters describing unparametrised geodesies is 6. Invoking the null condition, we see that 
the retraction projection of unparametrised null geodesies will, in general, constitute a 5-parameter 
family of curves in K 3 - in the special case where the conformal retraction is a static Killing vector, 
this is a 4-parameter family. 

On the other hand, let us consider a set of curves on some open set U C W 1 described by a 
system of third order ODEs. If we write this in an unparametrised way - as a set of (n — 1) third- 
order ODEs in terms of one of the coordinates - then we see that the integral curves of this system 
will lift to a foliation of the jet bundle J 2 (U, R) which is (3n — 2)-dimensional. 
Hence the unparametrised integral curves of a system of third-order ODEs in n dimensions consti- 
tutes a (3n — 3)-parameter family of paths. When n = 3, for example, we will have a 6-parameter 
family of such curves which is consistent with our results above. 

Hence, the set of projected null geodesies of the Kastor-Traschen metric, g, will constitute a 5- 
dimensional subset of the 6-dimensional family of unparametrised integral curves of ^ (except in 
the static case). 

So, a natural question arises: Given that we construct a third order system (such as ([8])) for which 
the projected null geodesies form a proper subset of the set of integral curves then to what do the 
other integral curves correspond? 

The system Q does not help us to answer this question but we can derive a different third order 
system which will. For convenience and clarity on this point, let us write our system of equations 
describing null geodesies of the Kastor-Traschen metric Q in three-dimensional vector notation 
i.e, 

x = (W - (W.x)x) ^(xx(Wx x)) . 



8 



Now let us consider a modification of this equation by adding an orthogonal term on the right-hand 
side, that is 

i= KTtf (ix(Wxx)) + A(ixW) (14) 

where A is a constant. Clearly, null geodesies satisfy this equation for A = 0. More interestingly, 
there is a six-parameter family of curves which satisfy this equation for some value of A. Hence, we 
might expect these curves to be the integral curves of some third order system in three dimensions 
which is independent of A. 

First from (14), we can derive the following equations by taking specific scalar products: 

2 



x.vv 

A 



V + cT 

2 

V + cT 
x.(x x W) 

|x x W| 2 



W.x + A(xx W).x 



VV x x| 



(15) 



Differentiating (14) and simplifying using the first two equations of (15), we derive the system of 

fcNV 



third order ODEs 



xx - - 



XV1 (xW)x-2c(xx(Vyxx))| X ' V1 



+ A 



2 |x x W| 2 
x.W 



X X 



x x W| 2 

(x.VV)x x VV + ((x x VV).x)x + x x VV + x x 



\ ds 



For A = 0, this system of equations reduces to M with the 



V+cT 



term replaced using &7n as ex- 



pected. We can eliminate A from this equation using the third equation of (15). Notice then that 
the vanishing of A coincides with x.(x x VV) = i.e, the vectors x, x and VV lie in the same plane. 
In particular, for the one centre case, the projections of null geodesies lie in a plane through the 
origin (centre). Overall, we have the following result 



Proposition 5.1 If c ^ 0, the retraction projection of the set of null curves satisfying (14) for 
some value of A coincides with the set of integral curves of 



3 x.W 



x = —xx— 



2 |x x VVl 



r (x.W)x - 2c(x x (VV X *)) 



x.W 

|x X W| 2 



fdVV 



x x 



\ ds 



x x 



x.(x x VV) 
|x x VV| 2 



x.W 

|x x W| 2 



(x.VV)x x W + ((x x W).x)x + — (x x VV) 

ds 



(16) 



Furthermore, the projected null geodesies are precisely the integral curves of this system for which 
x.(x x VV) = 0. 



Proof. As we have shown, any integral curve of (14) satisfies (16). To verify the reverse inclusion, 
we just need to consider the initial data unique to one integral curve 7 of ( 16 ) which will be given 
by seven parameters - three for initial position, two for initial unit velocity and two for initial 
acceleration (perpendicular to the velocity vector). By varying the values of T and A, it is clear 
that there is an integral curve of (14) with the same initial data and its projection necessarily 
coincides with 7. □ 

The proposition doesn't work for c = as we cannot use T as a parameter for the initial 
acceleration data in the above proof. 
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Magnetic Flow 



The addition of this extra A term may seem a little ad hoc here but is actually a sensible choice 
when we see the proof of this proposition i.e, we need to add a term orthogonal to x but not in 
the direction of x x (VV x x). Furthermore, this system of equations (14 1 can be interpreted as 
describing a magnetic flow in the background of the Kastor-Traschen metric with magnetic field 
B oc VV. This is precisely the magnetic field we saw in ^ when discussing the freedom in the 
2-form F and so these additional integral curves occupy a significant role in the geometry of the 
Kastor-Traschen metric. 



5.1 A Solution in the One-Centre Case 

Let ip — Amc and define a curve in the plane a; 3 = by 

x i (s)= f pa cob t^^log^s) j ,yssin f v/l ~ y ' 2 log(ys) J ,0j . (17) 



Then this curve satisfies hjki^x k = 1 and is an integral curve of the system of ODEs (16 1 for 
V = S. Furthermore, since it lies on a plane through the origin, we know, by Proposition 5.1, 
that it must be the retraction projection of a null geodesic of g. 

We can plot this curve in the plane and realise that it is just a reparametrised logarithmic spiral. 
This example is motivated by work in section 6 where we will show how to derive analytic expressions 



Figure 1: Logarithmic Spiral with ip = 0.1 



for the projected null geodesies of the one-centre K-T metric and thus, the integral curves of ( 10 1 



which lie on a plane through the origin. We should note here, that the limit c — > 0, for this example, 



is ill-defined. By L'Hopital's Rule, both expressions in (17) tend to zero, in this limit, for any value 
of s. 



6 Geodesies obtained from Extremal RNdS 

In this section, we will make use of a diffcomorphism between the one-centre Kastor-Traschen and 
the extremal Reissner-Nordstrom deSitter metrics. The advantage of this is that analytic solutions 
for the null geodesic equations of the RNdS metrics are well known [T3] . We will show how to derive 
these solutions which will enable us to obtain analytic solutions in the Kastor-Traschen coordinates 
and plot the retraction projection of the null geodesies in some cases. 
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6.1 RNdS Transformation 



We begin with the special Kastor-Traschen metric with potential V = ^ and now use spherical 
polar coordinates (|x| = R) to represent the flat metric h i.e, 

1 / 777 \ ^ 

9 = ~ 2 dT2 + \r +cT ) (dR 2 + R 2 (d9 2 + sin 2 9d<f> 2 )). (18) 

Assuming that c^Owe can make the coordinate transformation 

R = e~ cS , T= r -^e cS . (19) 
c 

If we choose t such that dt = dS + ^^^dr, where A r = (r — m) 2 — c 2 r 4 , then the metric becomes 

A r 2 

g = ^dt 2 + —dr 2 + r 2 (d9 2 + sin 2 dd<j> 2 ). (20) 

r z A r 

The analogous coordinate transformation in the Riemannian case was given in [5] . Here, the result- 
ing metric is the Reissner-Nordstrom deSitter spacetime in the extremal case with charge Q = rn 

A 



and c = ±J| where A is the cosmological constant. From now on, we will refer to t as static 



time in order to distinguish it from the time T. We also take c > as the definition of (201 is 
invariant with respect to c -> — c and we can compensate for it in the Kastor-Traschen case by 
sending T -> -T. 

Null and timelike geodesies of black hole spacetimes with cosmological constant were studied ex- 
tensively in [T3]. In particular, the author discussed the different types of orbits possible for 



Reissner-Nordstrom metrics and showed how to derive analytic formulae for the geodesies. The 

d_ 
dt 



RNdS metric admits the timelike static Killing vector Mr and we can plot the projection of the 



null geodesies to the space of orbits of Jj, analogous to what was done in [14], which we call the 
static projection. We then use this information to plot the null geodesies in the Kastor-Traschen 
metric projected along the conformal retraction 0, which, by Proposition 5.1, will be solutions of 



the system (16) and lie on a plane through the origin. 



Null geodesies of the RNdS metric ( 20 1 satisfy 



A r .„ r 2 



- -^-i 2 + —f 2 + r 2 (9 2 + sin 2 6 </> 2 ) = 0. (21) 
r z A r 

The Euler-Lagrange equation for 9 gives 

^-(2r 2 9) = 2r 2 sin6>cos6»</i> 2 , 
ds 

where s parameterizes the curves. By a choice of axes, we can set the initial conditions to be 
9 = 0, 9 = f, which results in motion in the equatorial plane - this coincides with the results of 
Proposition 5.1. Similarly, from the Euler-Lagrange equation for (f> and t we find that 

r 2 4> = $ 

2A r . 

T^t = -2E 
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where $ and E are constants. Let us focus our attention on non-radial geodesies and assume that 



$ > so that the geodesies are traced out in the direction of increasing (p. The null equation (21 1 
implies that 



or 



E 2 



-E 2 



([)- 







A, 



E 2 - V c 



eff- 



(22) 



Here V e // is the effective potential of the system which we can plot as a function of r. 




Figure 2: Plot of Veff as a function of r 



We note here that physically acceptable regions for null geodesies are those for which E 2 > V e ff. 
From the diagram, it is clear that we have two kinds of orbits - bound orbits (where r oscillates 
between two boundary values) and unbound flyby orbits (where r starts at oo, then approaches 
a periapsis and goes back to oo). This is consistent with the results of [2] and we use the same 
terminology. 



6.2 Circular Orbits 

As was observed in |15j . the extremum values of V e ff occur at r = m and r = 2m, independent 
of the cosmological constant (in this case, independent of c), resulting in circular orbits. The 
minimum of this function will not always be non-negative and, since E 2 > 0, it will not be attained 
for some values of c. Indeed, unless c = 0, we will not get a circular orbit at r = m. However, the 
transformation for R and T does not behave well in this limit and therefore does not allow us to 
see what this circular orbit corresponds to in the Kastor-Traschen coordinates. 
On the other hand, we will observe a circular orbit at r — 2m as long as 4mc < 1 (beyond this, the 
local maximum drops below the axis in Figure 2). For this solution, there will be some factors of c 
in the static time variable t and hence, when we make the coordinate transformation to R and T, 
the orbit in the retraction projection will be dependent on c. In fact, when we do this, the resulting 
orbit is just essentially a constant multiple of that given section 5 - a logarithmic spiral in some 
hyperplane of the space of orbits of the conformal retraction O, which passes through the origin. 



6.3 Horizon Structure and Nature of Orbits 

Clearly, for E 2 > V max , all orbits are unbound in the RNdS coordinates whereas for E 2 < V max , 
we can get both bound and unbound orbits. If m = then this graph becomes V e // = + const, 
and there are no bound orbits. 
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Using the chain rule and the E-L equation for <fr, we can rewrite equation ( 22 ) as 

2 



dr 
dcj> 



{k 2 + c 2 )r 4 - (r - to) 2 (23) 



where k = Also, 



dt _ r 4 
deb _ A r 



Note here that roots of A r (which we inevitably cross for some orbits) will cause infinities in ^ 
and in t itself - this will lead to a null geodesic tracing out a finite path in an infinite amount of 
static time (see analysis). Radii at which V e // = 0, or equally A r = 0, correspond to horizons. In 
particular, if: 

4mc < 1; there are three horizons (two black hole horizons and one cosmological horizon) 
with rbh- < fbh+ < Teh- The geometry is static for r < rbh-, rbh+ < r < r c h and corresponds 
to a black hole in a de-Sitter universe. 



4mc > 1; there is only one cosmological horizon. The geometry is static for r < r c h and 
corresponds to a naked singularity in a de-Sitter universe. 



The metric (20) has a singularity at r = which is covered by the Cauchy horizon r = rbh— in 
the case Amc < 1. The surface gravity at this inner horizon is larger than that at the cosmological 
horizon, in particular, 

- K ch = 8toc 3 > 

and so, by a result in [16] . the Cauchy horizon is unstable. Hence, some of the trajectories in the 
Reissner-Nordstrom-deSitter metric will be unphysical (in particular when r < rbh-)- Furthermore 
the case 4toc > 1 presents a possible violation of Penrose's cosmic censorship conjecture [T7] and 
may therefore also be unphysical. 

It was noted in [18] that, for the Reissner-Nordstrom deSitter metric in the extreme charge equal 
mass case, the Hawking temperature of the outer black hole horizon is the same, in magnitude, as 
that of the cosmological horizon endowing a notion of thermodynamic stability among all RNdS 
solutions. 

The plots of projected null geodesies can be obtained in the Kastor-Traschen framework by solving 



equation (23) and making a coordinate transformation (t, r) —> (T,R). As an example, we will 
perform this calculation for bound orbits in the three-horizon case (when 4toc < 1). The other 
curves can be similarly obtained by the reader but do not add much to the discussion. In the 
following example, we highlight when a trajectory is physical or when it is purely of mathematical 
interest. 

Bound Orbits in Kastor-Traschen with Three Horizons 

This case corresponds to the inequalities 

-1 + Vl +4mA 1 - \/l - 4toA 

4mA < 1 , < r < . 

2A ~ " 2A 



For bound orbits, any solution of (23) will oscillate between the two extremal values for r given 



above. Beginning with this equation, we can rearrange to get 



dv 

±1 , ^==0 + 7, (24) 
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where 7 is a constant of integration. With the given bounds on r, we can integrate the left hand 
side of this equation and obtain 



2F ^arcsin 



2Vl-4mA(-l-2rA+Vl+4roA) 



(l-2rA+v / l-4mA)(-2+v / l+4mA+ v / l-4mA) / ' 2 2 v / l-16m 2 A 2 



(1 - 16m 2 A 2 )V4 

where F is the elliptic integral of the first kind. We can solve this equation for r = r((f>) to get 
sin 2 ( J T 0))(1 + Vl-4mA)(Vl+4mA + yfl - 4mA - 2) + 2-y/l - 4mA(l - s/l + 4mA) 



'(</>) = 



where 



2A(sin 2 (J T ((/.))(Vl + 4mA + \/l - 4mA - 2) - 2-y/l - 4mA) 

1 



, (1 - Wm^ 2 ) 1 ^ ! 
J T (</>) = Jac =F — W + o - 



2 2-y/l - 16m 2 A 2 



and Jac is the Jacobi Amplitude for the elliptic integral (i.e, F(a,b) = c => a = Jac(c, b)). To 
obtain the oscillatory solution we use J(_i\n(0) whenever 



2F 



( "- 1 - 7 r I _ 1 

V 2 ^' 2 2\/l Ir T( 



(1 - 16m 2 A 2 ) 1 /4 



o p [ n 1 1 

16m 2 X 2 ) , ^ I 2 ' 2 2 x /I Ir T( 

> + 7 > -■ 



6m 2 A 2 



(1 - 16m 2 A 2 ) 1 /4 



We can now plot the static projection of the null geodesies of the RNdS metric. For this purpose, 



we take the values m = 1, k — g, c — | and 7 = 0. 





Figure 3: Static projection of null geodesic, < 
6 < 2tt 



Figure 4: < < IOtt 



As we mentioned before, it is important that we be careful here with respect to the range of the 
static time coordinate. In the above plots, the function r will, at several stages, cross a value for 
which A r = and satisfy r < r^-, where the trajectory is unphysical. This is reflected by the 
fact that each horizon crossing leads to an infinity in the static time t. For example, if we take 
a segment of this orbit which passes from the r^h- to rt,h+i we obtain the following plots for the 
geodesic itself and t as a function of (j> on this range. 
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r sin 




rcos^ 




Figure 5: Null geodesic traversing between black 
hole horizons and rf,h+- 



Figure 6: Static time t as a function of ■ 
interval. 



on this 



Using the transformation ( 19 1 together with the subsequent one for S 1 , we can plot the retraction 
projection of this null geodesic in the Kastor-Traschen coordinates. Similarly, we can determine 
the time T as a function of 4> and we get the following plots: 




Rcos0 



Figure 7: Retraction projection of null geodesic 
of the K-T metric 



Figure 8: Time T as a function of <j> on this inter- 
val. 



As R — > 0, the time T — > oo. However, we see that as R approaches the finite positive value, T 
also approaches a finite value. This suggests that we can extend our geodesic in the direction of 
decreasing T and continue the curve in the retraction projection. As one would expect, this can 
be done and the extension can be constructed by considering the part of the geodesic in the RNdS 
coordinates which begins at Thh- j decreases to the minimum value of r and then increases to rj>h— 
again. Of course, this curve lies completely within the Cauchy horizon and so, it is an unphysical 
extension but we can discuss it mathematically nonetheless. 

(Note: For this construction to provide the correct extension of the null geodesic, we must set the 
static time t to be decreasing on this interval). 

Hence, we can plot the full retraction projection of this null geodesic in the Kastor-Traschen coor- 
dinates. We include a plot of the time T here as a function of <p to demonstrate that it is, indeed, 
the projection of the whole geodesic - here the null geodesic begins at the spatial origin R = at 
T — » — oo, traces the curve in the direction of increasing and returns to the origin as T — > +oo: 
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Conformal Diagram 

To gain a better understanding of the geometry here, let us plot the trajectory of this null geodesic 
in the Penrose-Carter diagram of the spacetime. The authors in |19j have constructed this diagram 
and have highlighted the region covered by the "cosmological" (Kastor-Traschen) coordinates. We 
give a copy of this diagram and include a null geodesic which runs from a point at r < r^h- to the 
outer black hole horizon rbh+- 




Figure 11: Penrose-Carter diagram of extremal Reissner-Nordstrom deSitter spacetime. The region 
bounded by the red dotted line represents a single (T, R) chart. 



The authors of [19] make the point that the cosmological coordinates "smoothly cover the entire 
region from r — to r = oo." From the diagram, we see that a single (T, R) chart (bounded by the 
dotted red line) covers four (t, r) charts encompassing all three horizons. As the light ray crosses 
the unstable Cauchy horizon, the static time t — > oo. However in the Kastor-Traschen coordinates, 
we only have T — > oo as R — > and so, mathematically, the Kastor-Traschen metric does not see 
a time singularity at the inner Cauchy horizon. 
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7 Two-Centre Case 



Here, we make some remarks about the Kastor-Traschen metric with two black holes (N = 2). 
Analytic expressions for null geodesies are much more difficult to obtain but we can make some re- 
marks about the general theory. First of all, let us write our Kastor-Traschen metric g in cylindrical 
polar coordinates. 

! dT 2 + (V + cT) 2 (dp 2 +p 2 d(t) 2 + dz 2 ) (25) 



(V + cT) 2 

so that, without loss of generality, the singularities are placed on the z-axis, equidistant from the 
origin. That is 

y _ "*i rn 2 



{p 2 + (z - w) 2 ) 1 / 2 (p 2 + (z + w) 2 ) 1 / 2 

where m± and m 2 are the black hole masses and 2w is the distance between the centres. For this 
metric, we have the following result regarding the retraction projection of null geodesies. 

Proposition 7.1. // the retraction projection of a null geodesic has the property that its initial 
position and velocity lie in a plane passing through the two centres, then the entire projected null 
geodesic lies in that plane. 

Proof. Planes which pass through both centres are characterised by the condition <f> = constant. 
Then, the retraction projection of a null geodesic lies in such a plane if and only if (f) = at all 
points on the curve. Therefore, we can verify the proposition by showing that, at any point where 



:itv 

= 0, we have </> = 0. But, from equation rtsh , null geodesies of Kastor-Traschen satisfy 

2 . 2 fdV . dV \ ■ 

~P+ -7—Z (j) 



r p V + cT V<9p dz 

and the result follows. □ 

Now let us impose the additional condition mi = m 2 = M. Then, we discover another fixed 
plane of null geodesies. 

Proposition 7.2. If the retraction projection of a null geodesic has the property that its initial 
position and velocity lie in the plane passing orthogonally through the midpoint of the line segment 
joining the two centres, then the entire projected null geodesic lies in that plane. 

Proof. Clearly, the plane in question is given, in cylindrical polar coordinates, by z = and any 
null geodesic which lies completely in this plane will satisfy the condition i = at all points on 
the curve. Therefore, if the initial conditions z — 0, i = imply that z = initially, then the 
proposition is proved. Again, by equation ([5]), null geodesies of g satisfy 

dV fdV . dV 



-p + ~^r z 



and the result follows because 

dV 

dz 



V + cT \dz \dp^ dz 



Ma Ma _ 



(p 2 +w 2 ) 3 / 2 (p 2 +w 2 ) 3 / 2 

□ 
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7.1 Third Order System describing Null geodesies 



In section 5, we discovered that the third order system used to describe the retraction projection 
of null geodesies of the Kastor-Traschcn metric was not uniquely defined and by considering a new 
formulation we could obtain an interpretation of the entire set of integral curves. Here, let us take 
yet another system of third order ODEs, the integral curves of which contain the projected null 
geodesies of the one-centre Kastor-Traschen metric by eliminating the x l term from the system 
( |l0] ), using ([9) and Q to obtain 



Ixl 2 ^ 



2mc 



3(x.x) 



|x| 2 2(|x| 2 - (x.x) 2 ) 



(26) 



Then, we have the following result 



Proposition 7.3. Any integral curve of the system of ODEs (26), lies in a plane. Furthermore, 
if c =/= 0, such an integral curve will coincide with a projected null geodesic of the Kastor-Traschen 
metric g if and only if this plane passes through the origin. 



Proof. We construct the Frenet-Serret frame for a given integral curve of ( 26 1 



T = x 

Then, the Frenet-Serret formulas give 



N 



x 

W\ 



B = T x N. 



T = «N, 
N = -kT- 



tB, 



where k and r are the curvature and torsion of the curve, respectively. The first of these equations 



gives us k = |x|. Then, if we rewrite our system (26) in this frame, we obtain 

N = -kT 

and r = 0, necessarily. Hence, a given integral curve of this system must lie in a plane. 

We have already established that the projected null geodesies in the one-centre case will lie in 

a plane passing through the origin. To prove the "if" part of the proposition, we note that the 



initial data of an integral curve of (26) which lies on a plane through the origin is specified by 
six parameters - three for initial position, two for initial velocity (since it is unit in the arclength 
parametrisation) and one for the acceleration (in the plane of the position and velocity vectors, 
perpendicular to the velocity). Then, this curve is a projected null geodesic of the Kastor-Traschen 
metric if there exists a null geodesic with the same initial data. But, we can see that this is the 
case by analysing equation ([9]). Clearly, we can specify initial position and unit velocity vectors as 
we please. Then, the initial acceleration vector lies in the same plane perpendicular to the velocity 
and we can specify its magnitude by choosing the appropriate value of T. □ 



So, each integral curve of ( 26 ) lies in some plane but unlike the case of ( 16 1 , these integral curves 



have no obvious interpretation in terms of the Kastor-Traschen metric unless this plane passes 
through the origin. However, as we observe, they do arise as the projections of null geodesies for 
the two-centre case - in particular, those outlined by Proposition 7.2. 
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To see this let us consider the original expression for the Kastor-Traschen metric ([T]) with h as the 
flat metric in Cartesian coordinates with the potential written as follows 



rn 2 



V = 1 

|x — w| |x + w| 

where w = (w l ) is a fixed vector. For ease of notation, let us define the vectors 

™1 , i . ..is m 2 



X{ = 



:(x l + W % ). 



Then, null geodesies of the Kastor-Traschen metric are integral curves of the following system of 
third order ODEs 



3(X 1 +X 2 ).x(X 1 +X 2 ).x 
2 (Xi + X 2 ) . (Xi + X 2 - x) 

(Xi+X 2 ).x /3|x — w' 



-2c(Xi+^-((X 1 +X 2 ).x)x i 



(X x + X 2 ) . (Xi + X 2 - x) 



mi 



(Xlx)(i; - (XlxJx*) + 3|X + W| (X 2 .X)(^ - (X 2 .x)i<) 



TO 2 



This system is difficult to analyze, in general, but now let us assume that mi = m 2 = M and 
restrict attention to null geodesies which lie on the plane passing through the origin, orthogonal to 
the line between the two centres. 

Then |x + w| = |x — w| and since the acceleration and velocity vectors are perpendicular to w, we 
have w.x = w.x = 0. By making some simplifications using the geodesic equation Q, as in the 
one-centre case, we can replace the system above by 



x l = -\x\ 2 ± l 



Ulc 



- 3(x.x) 



2 x- 



(X.X)2) 



(27) 



Now we notice that this is precisely the system of third order ODEs ( 26 1 for the single centre case 



with black hole mass m = 2M, where |x — w| represents the distance from the centre. 
Hence, we have proved the following proposition: 

Proposition 7.4. Every null geodesic of the two-centre Kastor-Traschen metric which lies com- 
pletely in the plane passing orthogonally through the midpoint of the line segment joining the centres 



coincides with an integral curve of the system (26), with mass m = 2M , which lies in the plane a 



distance w from the origin, where 2w is the distance between the centres. 



Remark: This means that every integral curve of (26) can be realised as the retraction pro- 



jection of a null geodesic in either the one-centre or two-centre Kastor-Traschen metric, making all 
solutions physically relevant. 



7.2 Perturbation Analysis 

Now let us look at the stability of null geodesies in the z — plane by applying a small perturbation 
e <C a in the z-direction about the origin so that e and e are also small. Then substituting x 3 = z + e 



into (27) gives us the differential equation 



X / ■ "S n £ ~ 

1 1 I \z,z,z)— 



AMc 



- 3(x.x) 



2(|x- w| 2 - (x.x) 2 ) 



(z,z,z)=0 
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We can rewrite this as a coupled system of differential equations by choosing r\ = e and fi — e so 
that 



d 
ds 



I 4Mc 

|X-W|* 



x.x 

l*l a 



2(|x-w| 2 -(x.x) 2 ) 







(«,i,z)=0 



The stability of the system under small perturbations is determined by the eigenvalues of B. Given 
that the determinant of B is positive, (= |x| 2 ) we know that both eigenvalues have the same sign. 
If they are both positive then the system is unstable and if they are both negative then the system 
is stable. The mutual sign can be obtained from the trace of B and thus, we get the following result, 



Proposition 7.5: For the two-centre equal mass Kastor-Traschen metric, any geodesic which lies 
in the plane passing orthogonally through the midpoint of the line segment joining the two centres is 
stable under small perturbations normal to the plane at a point with given initial position, velocity 
and acceleration data if and only if 



- 3(x.x) 



4Mc / x.x 

|X- W| 3 ^|x]2 

at that point. Otherwise, it is unstable 



1 

|x - w| 



x.x 



2(|x- w| 2 - (x.x) 2 



< 



8 Unparametrised Projection of Null Geodesies in the One- 
Centre Kastor-Traschen Solution 



As final note, we will rewrite the system of third-order ODEs (26 1 in unparametrised form. By do- 
ing so, we get a purer notion of the set of projected null geodesies (free of parametrisation) and can 
make contact with the work in |20j where the author has explicitly derived differential invariants 
for systems of third order ODEs. 

To start with, let us relabel our coordinates x % = (z,x^) with (3 = 2,3. If we let ' denote differen- 
tiation with respect to z then we have 



■rP - fn-PY 



(x p )'z 



rP - (n.P\»& 



(x p )"z'' + (x^'z 



xP = (x^)'"i 3 + 3( X P)"zz + (xPyz. 



From the system (26), we can eliminate the "z term to obtain a pair of third order expressions 

2mc / x.x' 



(xPf'z 3 +3(x^)"iz 



3(x.x) 



1 



2(|x| 2 -(x.x) 2 ) 



(x )"z 2 . (28) 



This system can be simplified even further to eliminate the factors of z and z. First, let us use the 
following convention 



1 



o 



Then, using the arc-length parametrization condition hjkX J x k = 1, one can show that 



|u| 2 i 2 = 1,2 = 



u- 



(29) 
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Using this, it can also be shown that 



(u.a) : 



(x.u)(u.a) 



Then, we can use the expressions in ( 29 1 to eliminate z terms and subsequently powers of z in ( 28 1 



(This can be done by multiplying each term by the appropriate power of u z to give every term 
the same "weight" in terms of powers of z). The resulting expression will give us a pair of third 



order ODEs whose integral curves are the unparametrized curves of the system ( 26 1 



u.a 2mc|u| 3 / |u| 2 (x.a) — (x.u)(u.a) 
|u| 2 |a| 2 - (u.a) 2 



3(x.u) 



1 |u| 2 (x.a) - (x.u)(u.a) 
x| 2 2(|u| 4 |x| 2 - |u| 2 (x.u) 



Or") 



can be obtained on request (but are too long to be included with this work). 



In the limit c — > 0, the second and third terms on the right-hand side of (30) vanish coinciding 



with the unparametrised equation for conformal circles. I have calculated some of the Medvedev 
invariants [20 can for the system ( 30 1 and for conformal circles in Mathematica and the formulae 



9 Conclusions 

As was demonstrated in pQ, physical phenomena observed in spacetimes which admit a timelike 
vector field with a certain mathematical property can often be better understood by looking at the 
projection of null geodesies to the space of orbits of this vector field. This structure has been well 
documented in the static and stationary cases and in this paper, we discussed a specific example of 
a metric admitting a timelike conformal retraction which is also a solution of the Einstein-Maxwell 
equations and so, an important GR example. 

The third order system ^ arises naturally to describe the retraction projection of null geodesies 
of the Kastor-Traschen metric. We have, however, demonstrated that there is a freedom in the 
definition of this ODE system and that a more useful analysis is obtained by considering instead 



the system ( 16 1 where the totality of integral curves can be interpreted as the projection of null 
curves in the Kastor-Traschen metric describing a magnetic flow in the background magnetic field. 
This endows a physical relevance to this system and it would be interesting to probe its relevance 



further. Using this formulation, we've characterised those integral curves of (161 which coincide 
with the retraction projection of null geodesies. 

For the one-centre K-T solution, the projected null geodesies arc identified as those which lie on 



a plane through the origin. However, in this case, there is another third order system ( 26 ) whose 
integral curves arise as a distinguished subset of the projected null geodesies of the two-centre 
metric for some value of the distance between the centres, 2w, with masses mi = m% = M. This 
analysis of null geodesies appears to be a step further than has been seen thus far but extracting 
more analytic solutions is far from easy. 

There is a consistent physical intuition here if we consider what happens when w — > 0. We should 
expect to reproduce the retraction projection of some subset of the null geodesies for the one-centre 
Kastor-Traschen metric with black hole mass m = 2M and this is precisely what happens. In 
fact, we obtain all of the projected null geodesies because of the inherent spherical symmetry that 
accompanies this limit. 

One final point that we should stress here is that some of the physical properties of null geodesies 
obvious in the projection along one type of vector field can be obscured in the projection along 
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another. A clear example of this point can be seen in the one-centre Kastor-Traschen metric which 
we know, via the transformation to extremal RNdS coordinates, admits both a timelike conformal 
retraction and a timelike static Killing vector field. Light rays project down to unparametrised 
geodesies of the optical metric associated to the static Killing vector field and it can be shown that 
for different values of the cosmological constant A, the resulting optical metrics are projectively 
equivalent. One consequence of this is the fact that the differential equations governing light rays 
in these spacetimes are also independent of A. This phenomenon is not evident, however, when we 
consider the retraction projection of light rays in Kastor-Traschen coordinates. There is a clear c 
(or A) dependence in the system of ODEs ( 26 1 and this even carries over to the equations governing 
the unparamterised curves (30). 

Nonetheless, this is an interesting and physically relevant area of study and this paper encourages 
several open questions. It would be interesting to find analytic solutions to (16) in a more general 
case and to say something more concrete about the Kastor-Traschen metric with arbitrary V. 
Furthermore, we may ask about the properties of the null geodesic structure of an arbitrary metric 
which admits a timelike conformal retraction. 
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